Generalized form factors, generalized parton distributions 
and the spin contents of the nucleon 

M. Wakamatsu and Y. Nakakoji 1 

1 Department of Physics, Faculty of Science, 
Osaka University, 
Toyonaka, Osaka 560-0043, JAPAN 

Abstract 

With a special intention of clarifying the underlying spin contents of the nucleon, we investi- 
gate the generalized form factors of the nucleon, which are defined as the re-th x-moments of the 
generalized parton distribution functions, within the framework of the chiral quark soliton model. 
A particular emphasis is put on the pion mass dependence of final predictions, which we shall 
compare with the predictions of lattice QCD simulations carried out in the so-called heavy pion 
region around ~ (700 ~ 900) MeV. We find that some observables are very sensitive to the 
variation of the pion mass. It will be argued that the negligible importance of the quark orbital 
angular momentum indicated by the LHPC and QCDSF lattice collaborations might be true in 
the unrealistic heavy pion world, but it is not necessarily the case in our real world close to the 
chiral limit. 
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I. INTRODUCTION 



Q. 



The so-called "nucleoli spin crisis" raised by the European Muon Collaboration (EMC) 
measurement in 1988 is one of the most outstanding findings in the field of hadron physics 

0,0. 

The renaissance of the physics of high energy deep inelastic scatterings is greatly 
indebted to this epoch-making finding. Probably, one of the most outstanding progresses 
achieved recently in this field of physics is the discovery and the subsequent research of 
completely new observables called generalized parton distribution functions (GPD). It has 
been revealed that the GPDs, which can be measured through the so-called deeply- virtual 
Compton scatterings (DVCS) or the deeply- virtual meson productions (DVMP), contain 
surprisingly richer information than the standard parton distribution functions 

Roughly speaking, the GPDs are generalization of ordinary parton distributions and the 
elastic form factors of the nucleon. The GPDs in the most general form are functions of 
three kinematical variables : the average longitudinal momentum fraction x of the struck 
parton in the initial and final states, a skewdness parameter £ which measures the difference 
between two momentum fractions, and the four-momentum-transfer square t of the initial 
and final nucleon. In the forward limit t — > 0, some of the GPDs reduce to the usual quark, 
antiquark and gluon distributions. On the other hand, taking the n-th moment of the GPDs 
with respect to the variable x, one obtains the generalizations of the electromagnetic form 
factors of the nucleon, which are called the generalized form factors of the nucleon. The 
complex nature of the GPDs, i.e. the fact that they are functions of three variable, makes it 
quite difficult to grasp their full characteristics both experimentally and theoretically. From 
the theoretical viewpoint, it may be practical to begin studies with the two limiting cases. 
The one is the forward limit of zero momentum transfer. We have mentioned that, in this 
limit, some of the GPDs reduce to the ordinary parton distribution function depending on 
one variable x. However, it turns out that, even in this limit, there appear some completely 
new distribution functions, which cannot be accessed by the ordinary inclusive deep-inelastic 
scattering measurements. Very interestingly, it was shown by Ji that one of such distributions 
contains valuable information on the total angular momentum carried by the quark fields 
in the nucleon 

0-0. 

This information, combined with the known knowledge on the 
longitudinal quark polarization, makes it possible to determine the quark orbital angular 
momentum contribution to the total nucleon spin purely experimentally. 
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Another relatively-easy-to-handle quantities are the generalized form factors of the nu- 
cleon ^1,0], which are given as the non-forward nucleon matrix elements of the spin-n, 
twist-two quark and gluon operators. Since these latter quantities are given as the nu- 
cleon matrix elements of local operators, they can be objects of lattice QCD simulations. (It 
should be compared with parton distributions. The direct calculation of parton distributions 
is beyond the scope of lattice QCD simulations, since it needs to treat the nucleon matrix 
elements of quark bilinears, which are nonlocal in time.) In fact, two groups, the LHPC col- 
laboration and the QCDSF collaboration independently investigated the generalized form 
factors of the nucleon, and gave several interesting predictions, which can in principle be 
tested by the measurement of GPDs in the near future 0]. Although interesting, there 
is no a priori reason to believe that the predictions of these lattice simulations are realistic 
enough. The reason is mainly that the above mentioned lattice simulation were carried 
out in the heavy pion regime around m n ~ (700 ~ 900) MeV with neglect of the so-called 
disconnected diagrams. Our real world is rather close to the chiral limit with vanishing 
pion mass, and we know that, in this limit, the Goldstone pion plays very important roles 
in some intrinsic properties of the nucleon. The lattice simulation carried out in the heavy 
pion region is in danger of missing some important role of chiral dynamics. 

On the other hand, the chiral quark soliton model (CQSM) is an effective model of 
baryons, which maximally incorporates the chiral symmetry of QCD and its spontaneous 
breakdown 

Q,Q- (See Q Q 

for early reviews.) It was already applied to the physics 
of ordinary parton distribution functions with remarkable success For instance, 

an indispensable role of pion-like quark-antiquark correlation was shown to be essential 
to understand the famous NMC measurement, which revealed the dominance of the d- 



quark over the u-quark inside the proton |39j, 43], 44 1. Then, it would be interesting to 
see what predictions the CQSM would give for the quantities mentioned above. Now, the 
main purpose of the present study is to study the generalized form factors of the nucleon 
within the framework of the CQSM and compare its predictions with those of the lattice 
QCD simulations. Of our particular interest here is to see the change of final theoretical 
predictions against the variation of the pion mass. Such an analysis is expected to give some 
hints for judging the reliability of the lattice QCD predictions at the present level for the 
observables in question. 

The plan of the paper is as follows. In sect. II, we shall briefly explain how to introduce 
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the nonzero pion mass into the scheme of the CQSM with Pauli-Villars regularization. In 
sect. Ill, we derive the theoretical expressions for the generalized form factors of the nucleon. 
Sect. IV is devoted to the discussion of the results of the numerical calculations. Some 
concluding remarks are then given in sect.V. 



II. MODEL LAGRANGIAN WITH PION MASS TERM 



We start with the basic effective lagrangian of the chiral quark soliton model in the chiral 
limit given as 

C = ij(x)(i ft- MU^(x))ip(x), (1) 

with 



t/ 75 (x) 



3 i7sT-7T(a;)//x 



1 + 75 



U{x) 



75 



U\x). 



(2) 



2 v ' 2 

which describes the effective quark fields, with a dynamically generated mass M, strongly 



interacting with pions 



EL 



24( | . Since one of the main purposes of the present study is to see 



how the relevant observables depend on pion mass, we add to Cq an explicit chiral symmetry 
breaking term CJ given by [27 1 



C' = \flmltv f [U{x) + U\x)-2\. 



(3) 



Here the trace in (jSJ) is to be taken with respect to flavor indices. The total model lagrangian 
is therefore given by 

CCQM = £0 + C!. (4) 

Naturally, one could have taken an alternative choice in which the explicit chiral-symmetry- 
breaking effect is introduced in the form of current quark mass term as £' = — m ^V'- We did 
not do so, because we do not know any consistent regularization of such effective lagrangian 
with finite current quark mass within the framework of the Pauli-Villars subtraction scheme, 



as explained in Appendix of 
is given as 

with 



The effective action corresponding to the above lagrangian 
S eff [U] = S F [U) + S M [U], (5) 

S F [U] = -i N c Sp \n(i @ - MW 5 ) , (6) 
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and 

Sm[U] = J\fl ml try [U{x) + U\x) - 2] d'x. (7) 

Here SpO = / rf 4 xtr 7 try (x|0|x) with tr 7 and try representing the trace of the Dirac gamma 
matrices and the flavors (isospins), respectively. The fermion (quark) part of the above action 
contains ultra-violet divergences. To remove these divergences, we must introduce physical 
cutoffs. For the purpose of regularization, here we use the Pauli-Villars subtraction scheme. 
As explained in |22(, we must eliminate not only the logarithmic divergence contained in 
Sf[U] but also the quadratic and logarithmic divergence contained in the equation of motion 
shown below. To get rid of all these troublesome divergence, we need at least two subtraction 
terms. The regularized action is thus defined as 

s r e e / f [u} = s r F e9 [u} + s M [u], (8) 

where 

s r ; 9 [u] = s F [u]-j:c l sp[u]. (9) 

1=1 

Here Sp is obtained from S^t/] with M replaced by the Pauli-Villars regulator mass Aj. 
These parameters are fixed as follows. First, the quadratic and logarithmic divergence 
contained in the equation of motion (or in the expression of the vacuum quark condensate) 
can, respectively, removed if the subtraction constants satisfy the following two conditions : 



M 2 -J>A 2 = 0, (10) 

1=1 

m 4 -E^ = o. (11) 

i=i 

(We recall that the condition which removes the logarithmic divergence in Sf[U] just coin- 
cides with the 1st of the above conditions.) By solving the above equations for c\ and C2, 
we obtain 

which constrains the values of c\ and C2, once Ai and A2 are given. For determining Ai and 
A 2 , we use two conditions 

2 2 / A . \ 2 / A . \ 2 



N C M 

~Jtx 2 frf \.\f I "' \ \[ 
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and 

which amounts to reproducing the correct normalization of the pion kinetic term in the 
effective meson lagrangian and also the empirical value of the vacuum quark condensate. To 
derive soliton equation of motion, we must first write down a regularized expression for the 
static soliton energy. Under the hedgehog ansatz tz{x) = f n rF(r) for the background pion 
fields, it is obtained in the form : 

E r Zi C [F(r)} = E r p 9 [F{r)} + E M [F(r)), (16) 

where the meson (pion) part is given by 

E M [F{r)) = -flml J d 3 x [cosF(r) - 1], (17) 

while the fermion (quark) part is given as 

E r F e9 [F(r)) = E val + E:;°, (18) 



with 



E val = N C E , (19) 



N 



K P 9 = N c £ (E n - E^) -E^JVeE (Et - E^). (20) 

n<0 i=l n<0 

Here E n are the quark single-particle energies, given as the eigenvalues of the static Dirac 
Hamiltonian in the background pion fields : 

H\n) = E n \n), (21) 

with 

a ■ V 

H = — — + /3M [cos F(r) + ij 5 T ■ r sin F(r)}, (22) 
i 

while E^ denote energy eigenvalues of the vacuum Hamiltonian given by (|2*2*|) with F(r) = 
(or U = 1). The particular state |n = 0), which is a discrete bound-state orbital coming 
from the upper Dirac continuum under the influence of the hedgehog mean field, is called the 
valence level. The symbol J2 n <o i n (BOj) denotes the summation over all the negative energy 
eigenstates of H, i.e. the negative energy Dirac continuum. The soliton equation of motion 
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is obtained from the stationary condition of E r s l^ tic [F {r)\ with respect to the variation of the 
profile function F(r) : 

= Ji-E^Fir)] 

= 4vrr 2 {-M[S(r) sin F(r) - P(r) cos F(r)\ + f*m 2 w sin F(r) } , (23) 

which gives 

F(r) = arctan ( — . — — ) . (24) 
Here 5(r) and P(r) are regularized scalar and pseudoscalar quark densities given as 

S(r) = S val (r) + E S n (r) "E^E S^(r), (25) 

n<0 i=l iW n<0 

P(r) = P val (r) + J2 Pn(r) "E^E *£(r), (26) 



with 



^/d 3 * Ha,) 
47r 7 r 



S„(r) = d 3 x (n\x) 7 ° (aj|n>, (27) 



Pn(r) = -^J<?x(n\x)i'f'ftT'T ll ^ 2 rJ (x|n), (28) 

while (r) and P Ai (r) are the corresponding densities evaluated with the regulator mass 
A, instead of the dynamical quark mass M. We also note that S va i(r) = 5 n= o(r) and 
P va i(r) = P n= o(r). As usual, a self-consistent soliton solution is obtained as follows with 
use of Kahana and Ripka's discretized momentum basis j^jiHl- First by assuming an 
appropriate (though arbitrary) soliton profile F(r), the eigenvalue problem of the Dirac 
Hamiltonian is solved. Using the resultant eigenfunctions and their associated eigenenergies, 
one can calculate the regularized scalar and pseudoscalar densities S(r) and P{r). With use 
of these S(r) and P(r), Eq.(|24j) can then be used to obtain a new soliton profile F(r). The 
whole procedure above is repeated with this new profile F(r) until the self-consistency is 
attained. 



III. GENERALIZED FORM FACTORS IN THE CQSM 



Since the generalized form factors of the nucleon are given as moments of generalized 
parton distributions (GPDs), it is convenient to start with the theoretical expressions of 
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he unpolarized GPDs H(x,£,t) and E(x,£,t) within the CQSM. Following the notation in 
28[,|29|, we introduce the quantities 

,(1=0) _ fd\ iXx , , -r ( \n\ , , ( \n 



and 



M '' 7 I 2i e <P' S 'I* T n*[ T \\p,'), (29) 



M«f> 3 / £ <f (p'_ S '| 4- r 3 M (£j |p, .>. (30) 

Here, the isoscalar and isovector combinations respectively correspond to the sum and the 
difference of the quark flavors u and d. The relation between these quantities and the gener- 
alized parton distribution functions H(x,£,t) and E(x,^,t) are obtained most conveniently 
in the so-called Breit frame. They are given by 

3kl A k, 

M^ 0) = 25 s , s H < i="\x } i J t) - (a l ) s/s E^=°\x,^t), (31) 



3kl A k 

M$f> = 25,.H% = *\x,t,t) - IL—^^ix^t)- (32) 

M N 



where 



Ht 0/1) (x,U) - H^\x,U) + ^E^\x,U), (33) 
E^r m \x^,t) H^\x,i,t) + E^{x^,t). (34) 

These two independent combinations of H(x,£,t) and E(x,£,t) can be extracted through 
the spin projection of M.^ I= °/^ as 

H%\x,£,t) = \tr{M^}, (35) 
E$(x,Z,t) = iMN ^T Ab H- m M^}, (36) 

where "tr" denotes the trace over spin indices, while = A 2 — (A 3 ) 2 = —t — (— 2Mn£) 2 . 
Now, within the CQSM, it is possible to evaluate the right-hand side (rhs) of and (|3*Sj) . 



28]— 



30|, we do not repeat 



Since the answers are already given in several previous papers 
the derivation. Here we comment only on the following general structure of the theoretical 
expressions for relevant observables in the CQSM. The leading contribution just corresponds 
to the mean field prediction, which is independent of the collective rotational velocity Q of 
the hedgehog soliton. The next-to-leading order term takes account of the linear response 
of the internal quark motion to the rotational motion as an external perturbation, and 



consequently it is proportional to Q. It is known that the leading-order term contributes 
to the isoscalar combination of HE{x,£,t) and to the isovector combination of EM{x,£,,t), 
while the isoscalar part of //^(x, and the isovector part of EM(x,£,t) survived only at 
the next-to- leading order of Q (or of 1/N C ). The leading-order GPDs are then given as 

dz° 

n<0 

X (1 + 7° 7 3 ) e -*(*Wa e iA ' x e -*(*°/ 2 )ft> <& n (x) 



H { t°\x, 



(1=1), 

M 1 



2iM 2 N N c rdz 



aj 



x (1 + 7° 7 3 ) (r x A) 3 e" i(z ° /2)?53 e iA a; e~ i{z ° /2) ^ $ n i 



(37) 



(38) 



Here the symbol E n < denotes the summation over the occupied (the valence plus negative- 
energy Dirac sea) quark orbitals in the hedgehog mean field. On the other hand, the theo- 
retical expressions for the isovector part of He and the isoscalar part of Em, which survive 
at the next-to-leading order, are a little more complicated. They are given as double sums 
over the single quark orbitals as 



H 



(1=1), 



M N N r [df_ 
27 



12/ 



E 

v m=all,n<0(E m ^E n ) 



-iE n z' 







E 

n=all,m<0(E m y£E n ) 



-iE m z" 



E - E 



Id . 

+ Mat dx ^ 



-iE n z° 



ixM N z° 



m=all,n<0(E rn ^ 1 E n ) 

x (n|r a |m)(m|r a (l + 7°7 3 )e~ i(z ° /2)?53 e iA - :r e- i(20/2)?53 |n). (39) 



and 



. M 2 N N C rdf_ 
1 21 J 2tt 



E 

m=all,n<0(E m ^E n ) 



-iE n z 







E 

n=all,m<0(E m y£E n ) 



-iE m z 







E — E 

J-'n J -'m 



Id 



-iE n 



M N dx m=aUyn < ( Em ^E n ) 



JxM N z° 



x (n\r b \m)(m\ (1 + 7 7 3 ) e~ l{z ° ^ 



p3o6^q 



iA-a; „-i( z °/2)p 3 



n). (40) 



These four expressions for the unpolarized GPDs, i.e. (|S7|) ~ (|4H|l. are the basic starting 
equations for our present study of the generalized form factors of the nucleon within the 



CQSM. There are infinite tower of generalized form factors, which are defined as the n- 
th moments of GPDs. In the present study, we confine ourselves to the 1st and the 2nd 
moments, which respectively corresponds to the standard electromagnetic form factors of 
the nucleon and the so-called gravitational form factors. We are especially interested in the 
second one, since they are believed to contain valuable information on the spin contents of 
the nucleon through Ji's angular momentum sum rule j^Q]. For each isospin channel, the 
1st and the 2nd moments of H^\x,0,t) define the Sachs-electric and gravito-electric form 
factors as 

Gg,T 0/1) (t) = £ H«= 0/1) (x,0,t)dx, (41) 

and 

G^=o 0/1) (0 = £ x H { z =0/1) (x, 0, t) dx. (42) 

On the other hand, the 1st and the 2nd moments of E^/J (x,0,t) respectively define the 
Sachs-magnetic and gravito-magnetic form factors as 

Gfi*™ (t) = £ E^= 0/l) (x,0,t)dx, (43) 

and 

G^T (t) = £ x I?2 =0/1) (x, 0, t) dx. (44) 

In the following, we shall explain how we can calculate the generalized form factors based 
on the theoretical expressions of corresponding GPDs, by taking G^To (0 an< ^ ^e^o\^) as 
examples. Setting £ = and integrating over z° in (13 7|) . we obtain 

#i J=o) (a;,0,t) = M N N C J d 3 xe iA± x 

x J2 & n (x)(l + -yV)5(xM N -E n -p 3 )$ n (x). (45) 

n<0 

Putting this expression into (|41jh we have 

G%${t) = N c f d 3 xe iA ^ x +7°7 3 )$n(x). (46) 

It is easy to see that, using the generalized spherical symmetry of the hedgehog configuration, 
the term containing the factor 7°7 3 identically vanishes, so that G^^it) is reduced to a 
simple form as follows : 

G^7 0) (t) = N c f d 3 x £ *t (*) j (A ± x) Q n (x). (47) 
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Aside from the factor N c (= 3), this is nothing but the known expression for the isoscalar 
Sachs-electric form factor of the nucleon in the CQSM j^J]. 

A less trivial example is Gg~^ (t), which is defined as the 2nd moment of H^~°\x, 0, t). 
Inserting (J45j) into (|41|) and carrying out the integration over x, we obtain 



G%%\t) = ^N c f d 3 xe iA -- x £ *i(*)(l + 7V)(i^ + ft) *»(*)■ (48) 

M N n <0 

Using the partial-wave expansion of e iA± - x , this can be written as 
G%£\t) = -±-N c fd 3 xY: 47ri^(A ± ) 

Mm J , 



N J l, m 

x ^2^ n (x)j l (A ± x)Y lm (x){l + a 3 )(E n +p 3 )^ n {x). (49) 

n<0 



This can further be divided into four pieces as 



g£o 0) W = Eg 8) (50) 



i=l 



where 



with 



G% = W Nc I ^ ^ 4 ^m(A±)M;, (51) 



M x = J2 ®l( x ) Ji( A ± x ) Y im(x) E n ® n (x), (52) 

n.<0 

M 2 = & n {x)]i(A ± x)Y lm (x)a 3 E n <5> n (x), (53) 

n<0 

M 3 = ]T ^ n (x)ji(A ± x)Y lm (x)P3<l>n(x), (54) 

n<0 

Mi = E ^i(sc)ji(A±a;)yi m (x)a3P3^n(aj). (55) 

n.<0 

To proceed further, we first notice that, by using the generalized spherical symmetry, Mi 
survives only when I — m — 0, i.e. 

Mi ~ £ ^(x)jo(A ± x)^^E n $ n ^), (56) 

n<0 V47T 

which leads to the result : 

N r 

G i = XT / ^ £ ®n( x ) 3o(A±x) E n <& n (x). (57) 

M N J „<o 
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To evaluate G 2 , we first note that 

Y lm {x) xa 3 = E {lmlO\Xm) \Y t (x) x a]< A >. (58) 

A 

Here, the generalized spherical symmetry dictates that A must be zero, so that the rhs of 
the above equation is effectively reduced to 

- ^u^o^i^xq]' ). (59) 

This then gives 

G2 = ^ N cJ d 3 x4mY* (A ± ) (-^ £ < («)ji(A±rr) [Yi(x) x a]^ * B (*).(60) 
Owing to the identity 

Yo(A ± ) = JTp^onl) = 0, (61) 

we therefore find that 

G 2 = 0. (62) 
Next we investigate the third term G 3 . Using 

Y lm (x) x p 3 = E (/mlO|Am> [Y^x) xp]W ^= 5 U 5 m , [Fi(x) x p]<°), (63) 

a V3 

we obtain 

Gs = J^ N cjd 3 x4niY* (A ± ) (- E ®n( x )h{A±x) \Yi(x) x p] (0) (64) 

This term vanishes by the same reason as G2 does. The last term G4 is a little more 
complicated. We first notice that 

Y lm {x) a 3 p 3 = Y lm (x) E (1010|A0) [a x pf ] 

A 

~ E (1010|A0) (ZmA0|00) [Yi(x) x [a x p]W](°) 

A 

= 5 mfl (1010|/0) (/0/0|00) \Yi(x) x [a x p] w ] (0) , (65) 

which dictates that I must be or 2. Inserting the above expression into 1)5 5 j) . and using the 
explicit values of Clebsch-Gordan coefficients, G4 becomes 
1 Nc_ 



G, = --= ■ -f- fd 3 x £ $t .(«) jo(Axx) [a x p](°) $ n («) 
\/3 Mv •/ n < 

I rf 3 x E $t ( x )j 2 (A ± x) [Y 2 (x) x [a x pp)]C°> $„(*). (66) 



V47T iV c 



n<0 
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Using the identities 



[axp]<°> = -^"-P, (67) 
[y 2 (x) x [a x p]( 2 )](°> = [[Y 2 (£) x p]« x a](°), (68) 



G> = \- N " 



G 4 can also be written as 

o 1 , / rf3;r E ^fa) io(A ± x) a • p $ n (x) 
3 M N J n < 

+ • XT / d " x E j^s) [T 2 (x) x p]« x a]W (69) 

v 6 M;v J n < 

Collecting the answers for Gi, G 2 , G 3 and G 4 , we finally obtain 

Cg 20 0) (0 = -vf- / d 3 * £ *U*)io(A ± x)£ n <t> n (*) 
M Jv ^ n<0 



1 iV c 



3 Mjy 



V4vr A" c 
+ ' 



/* d 3 x ®l( x ) Jo(A±x) a • p $ n (au) 

^ n<0 

f / rf 3 X £ M*±X) M&) >< P] W >< «] (0) $ n(*)- (70) 

V6 Mat J n < 

Up to now, we have obtained the theoretical expressions for the isoscalar combination 
of the generalized form factors G^J^t) and G^^\t). For notational convenience, we 
summarize these results in a little more compact forms as follows : 

= r ^ =0) (*,0,t)rf* = ^E (n| Jo (A ± r)|n), (71) 
J - 1 „<o 



and 



1 



N c ]T E n {n\j (A ± r)\n) 



1 



+ N C J2 (n\j (A ± r)-a-p 

n<0 6 



n 



+ E Hj 2 (A ± r) [[y 2 (f) x p]W x a]<°> |n> . (72) 

v6 n < J 

As pointed out before, G^~^{t) is N c (= 3) times the isoscalar combination of standard 
Sachs-electric form factor of the nucleon. Analogously, we may call G^^\t) the gravito- 
electric form factor of the nucleon (its quark part), since it is related to the nonforward 
nucleon matrix elements of the quark part of the QCD energy momentum tensor. 
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The other generalized form factors can be obtained in a similar way. The isovector part 
of the generalized electric form factors survive only at the next-to-leading order of Q. They 
are given as 



1 /N c \ 1 

' p (m\\T\\n)(m\\j {A ± r)T\\n), (73) 



and 



31 V 2 / m>0 n < -Em -En 

= £ xHV =1 \x,0,t)dx 



M N 31 V 2 ) m>0 n < E m - £ n 
— - (m||j (Aj.r)T||n) 



+ (m||j (Aj.r) - (a -p)r||n) 



+ l^.{ m \\j 2 (A ± r)[[Y 2 (f) x p]« x a](°>T||n>J. (74) 

The isoscalar combination of the generalized magnetic form factors also survive only at the 
next-to- leading order of Q, so that they are given as double sums over the single-quark 
orbitals in the hedgehog mean field as 

M " /N f) E ^yHHWHI^frxaHW, (75) 



and 



G$$(t) = j\E£-°\x,0,t)dx 



1 fN r 



— I — -) — — (m||r||n) 

T \ 2 J ^ E — E 



m>0,n<0 ^ m 

j E m + E n ji(A ± r) ji(A±r) ,, ,1 , * 

x < (ml (r x ct) \ In) + (mi — L n) > . (76) 

^ 2 Aj_r ' A^r ' J 

We recall that (i) just coincides with the known expression of the isoscalar Sachs- 

magnetic form factor of the nucleon in the CQSM 45]. On the other hand, G^~ 2 °o(t) is 
sometimes called the gravitomagnetic form factor of the nucleon (its isoscalar part), which 
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we can evaluate within the QCSM based on the above theoretical expression. Finally, the 
leading-order contribution to the isovector part of the generalized magnetic form factors are 
given as 

G$*j(t) = I" E$ =1 \x,0,t)dx = -^f ■ N c J2(n\\ J -^^r ■ (r x a )\\n), (77) 

6 n<0 l\±T 



and 



(t) = J xE ( M =1 \x,0 } t)dx 



= -i-^i;{^(n||^r.(rxa)|| n ) + HI . L \\„) } .(78) 

Especially interesting to us are the values of the generalized form factors in the forward 
limit t — *■ 0. The consideration of this limit is also useful for verifying consistency of our 
theoretical analyses, since it leads to fundamental sum rules discussed below. We first 
consider the forward limit of G^7o (0- From (|71|). we find that 

Gtioit = 0) = C H { t°\x,0,0)dx = N c J2 I- (79) 

J ~ X n<0 

Subtracting the corresponding vacuum contribution, this reduces to N c (= 3). If we remem- 
ber the relation 

J H { E =O \x,0,0)dx = ^ H u+d {x, 0, 0)dx = J 1 J u+d {x)dx = N u + N\ (80) 

the forward limit of (J7T|) just leads to the sum rule : 

{t = 0) = N u + N d = 3, (81) 

which denotes that the sum of the u-quark and (i-quark numbers in the proton is three. 
Next we turn to the forward limit of G^~^\t), which gives 

G^it = 0) = ^- ( N c £ E n + \ N c (n\a ■ p\n) 1 . (82) 

[ n < 3 n < J 

It is easy to see that, after regularization and vacuum subtraction, the first term of the rhs 
of the above equation reduces to the fermion (quark) part of the soliton energy, i.e. E 7 p 9 
in (JIBJ). It was proved in that, in the CQSM with vanishing pion mass, the following 
identity holds : 

(n\a -p\n) = (83) 



n<0 
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In the case of finite pion mass, which we are handling, this identity does not hold. Instead, 
we can prove (see Appendix) that 

\N c Y,{n\cx-p\n) = E M . (84) 

That is, the second term in the parenthesis of rhs of eq.([82|) just coincides with the pion 
part of the soliton energy (or mass). Since the sum of the quark and pion part give the total 
soliton mass M N , we then find that 

G { E^\t = 0) = -L.M N = l. (85) 
M N 

In consideration of eq.(J72J), this relation can also be expressed as 

J^xH { z =0 \x,0,0)dx = J xH u+d (x,0,0)dx = J xf u+d (x) dx = (x) u+d = 1, (86) 

which means that the total momentum fraction carried by quark fields (the u- and <i-quarks) 
is just unity. This is an expected result, since the CQSM contains quark fields only (note 
that the pion is not an independent field of quarks), so that the total nucleon momentum 
should be saturated by the quark fields alone. 

Taking the forward limit of G^l\t), we are again led to a trivial sum rule, constrained 
by the conservation low. In fact, we have 

G^(t = °) = 7 (t) £ H|T||n)2 = 7 ■ 1 = x > (87) 

1 V 6 J m>0,n<0 Urn-tin I 

thereby leading to 

I^H { E =l \x,0,0)dx = ^ H u ' d {x, 0, 0) dx = J 1 J u - d (x)dx = N u -N d = l, (88) 

which denotes that the difference of the u-quark and the d-quark numbers in the proton is 
just unity. On the other hand, the forward limit of {t) leads to the first nontrivial sum 

rule as 



u—d 



G^- 1) (t = 0) = J 1 xH ( z- 1) (x,0,0)dx = ^ xf u ~ d (x)dx = (x) 
1 1 fN c \ _ 1 , „ „ v 

= — — — I — I \ m \\ T \\ n ) 

m n 1 V 6 J m>0jn < E m - E n 

x | m ~^ — -(m||r||n) + (m|| — (at ■ p) t \\n)\ . (89) 
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Since this quantity, which represents the difference of momentum fraction carried by the u- 
quark and the ci-quark in the proton, is not constrained by any conservation law, its actual 
value can be estimated only numerically. 

Next we turn to the discussion of the forward limit of the generalized magnetic form 
factors. First, the forward limit of Crj^io (0 gi yes 

Gfi§(« = 0) = -^f(^r) E ^p-HHInXmllrxalln), (90) 

1 V b J m>0,n<0 h ra ~ E n 

which reproduces the known expression of the isoscalar magnetic moment of the nucleon in 



the CQSM 



3- On 



the other hand, the forward limit of G^f^(t) gi yes 



<#> = <» = £ r^HHW 



m>0,n<0 m 

( E m + E n 

X 



| ^ m ^~ ^ n (m\ \r x at\ |n) + (ra||JD||n)| . (91) 
It was shown in j3] that the rhs of the above equation is just unity, i.e. 

G$${t = 0) = 1- (92) 
In consideration of (jSj), this identity can be recast into a little different form as 

i = y xE^ =0) (x,o,o)rfx 

= J^xE^ d (x,0,0)dx = J 1 x [H u+d (x, 0, 0) + £ u+d (x, 0, 0)] dx. (93) 
Assuming the familiar angular momentum sum rule due to Ji 

- J x[H u+d (x, 0, 0) + E u+d (x, 0, 0)] dx = J u+d , (94) 
the above identity claims that 

J u+d = \ (95) 

which means that the nucleon spin is saturated by the quark fields alone. This is again a 
reasonable result, because the CQSM is an effective quark model which contains no explicit 
gluon fields. The derived identity ()92j) has still another interpretation. Remembering the 
fact that G^2Q (t) consists of two parts as 

G { L:$(t)=A% d (t) + B£ d (t), (96) 
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En.(|52jl dictates that 

A« +d (0) + i? 2 V d (0) = l. (97) 
Since it also holds that (the momentum sum rule) 

G^S(0) = A^ d (0) = l, (98) 

it immediately follows that 

5 2 V d (0) = 0, (99) 

which is interpreted as showing the absence of the net quark contribution to the anomalous 
gravitomagnetic moment of the nucleon. 

Finally, we investigate the forward limit of the isovector combination of the generalized 
magnetic form factors. From eq. (|77|). we get 

G { L:${t = 0) = -^N c J2(n\\r-(rx a)\\n), (100) 

y n<0 

which reproduces the known expression of the isovector magnetic form factor of the nucleon 
in the CQSM. On the other hand, letting t — > in (|78|). we have 

siftS (* = o) = -^cE{^HI T -( rxa )ll n > + H|T-^||n)} 

J n<0 

1 a;E^ =1) (x,0,0)cix. (101) 

As shown in 29j. this sum rule can be recast into the form : 

l - J xE { ^ =1 \x, 0, 0) alx = J {1=1 \ (102) 

where j( I=1 ^ consists of two parts as 

jC=i) = Jf =v > + 5JV =1 \ (103) 

Here, the first part is given as a proton matrix element of the free field expression for the 
isovector total angular momentum operator of quark fields as 

Jf =1) = (p]\jf =1) \p]), (104) 



with 



(x x p) 3 + ^E 3 



L^ + lE^i). (105) 
18 



On the other hand, the second term is given as 



gj(i=i) = -M(— ) V(n|rsinF(r)7 [S-fT-f-S-r]|r2). 

18 ) „<0 



(106) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



The model in the chiral limit contains two parameters, the weak pion decay constant f n 
and the dynamical quark mass M. As usual, /„- is fixed to its physical value, i.e. f n = 
93 MeV. For the mass parameter M, there is some argument based on the instanton liquid 
picture of the QCD vacuum that it is not extremely far from 350 MeV |23j. The previous 
phenomenological analysis of various static baryon observables based on this model prefer a 
slightly larger value of M between 350 MeV and 425 MeV 3] 2^|. In the present analysis, 
we use the value M = 400 MeV. With this value of M = 400 MeV, we prepare self- 
consistent soliton solutions for seven values of m^, i.e. = 0,100,200,300,400,500, and 
600 MeV, in order to see the pion mass dependence of the generalized form factors etc. 
Favorable physical predictions of the model will be obtained by using the value of M = 
400 MeV and = 100 MeV, since this set gives a self-consistent soliton solution close to 
the phenomenologically successful one obtained with M = 375 MeV and m n = MeV in 
the single-subtraction Pauli-Villars regularization scheme 

We first show in FigQJthe soliton profile functions F(r) obtained with several values of 
m^, i.e. m n = 0,200,400, and 600 MeV. One sees that the spatial size of the soliton profile 
becomes more and more compact as the pion mass increases. 

We are now ready to show the theoretical predictions of the CQSM for the generalized 
form factors. Since the corresponding lattice predictions are given for the generalized form 
factors A^ d (Q 2 ) and B^ d (Q 2 ), which are the generalization of the standard Dirac and 
Pauli form factors, we first write down the relations between these form factors and the 
generalized Sachs-type factors, which we have calculated in the CQSM. They are given by 



A^ d it) 
"CO 

A u — d I 
/ii 



AU+d 

^20 



(t) 
■»o (*) 



L 10 

j^u—d 



L. 10 ft) + T 



G (i=o) 
G (h in 



G 



(i=l), 



t(/=o), 

r M,10 1 

t(/=0), 
r A/,20 > 



(*) /(1 + r), 



(t) /(1 + r), 



<%*\t)+TG%£>(t) /(1 + r), 



(107) 
(108) 
(109) 
(110) 
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m„ = MeV 
m„ = 200 MeV 
m„ = 400 MeV 
ITL = 600 MeV 



8.0 



10.0 



FIG. 1: The self-consistent soliton profile functions obtained with M = 400 MeV and m T = 
0,200,400, and 600 MeV. 



and 



sT(t) = 


- r (i=o) (f \ 




+ 


(111) 


B£*(t) = 




~ ^£,20^(0 




(112) 


B£*(t) = 




_ ^£,10^) 




(113) 


B u 2 , d {t) = 




~ ^£,20^) 


+ 


(114) 



where r = —t/AM%. We recall that A w (t) and B w (t) are nothing but the standard Dirac 
and Pauli form factors of the nucleon : 



A^ d (t) = Fl I=0 \t)=F?(t) + F?(t), 



B^ d {t) 



Ft°\t)=F 2 p (t) + F2 a (t), 
Fl I=1 \t) = Ff(t) - F?(t), 



B u w d {t) = Fl'- 1 \t)=F?(t)-F?(t). 



(115) 
(116) 
(117) 
(118) 



Since the lattice simulations by the LHPC and QCDSF collaborations were carried out 
in the heavy pion region around m n ~ (700 ~ 900) MeV and since the simulation in the 
small pion mass region is hard to perform, we think it interesting to investigate the pion 
mass dependence of the generalized form factors within the framework of the CQSM. For 
simplicity, we shall show the pion mass dependence of the generalized form factors at the 
zero momentum transfer only. We think it enough for our purpose because the generalized 
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form factors at the zero momentum transfer contain the most important information for 
clarifying the underlying spin structure of the nucleon. At zero momentum transfer, the 

relations between the generalized Dirac and Pauli form factors and the generalized Sachs- 
type form factors are simplified to become 

^io +d (0) = ^(O), (119) 

A^ d (0) = Gfi 0) (0), (120) 

Atf(0) = Gg^(0), (121) 

4»*(0) = Ggjfto), (122) 

and 

B u 10 +d (0) = G$$(0) - Gg 0) (0), (123) 

B^ +d (0) = G$$(0) - Gg o 0) (0), (124) 

*iV(0) = ^(0) - Gg^O), (125) 

B£*{0) = GS3?(0) - Gg'lO). (126) 
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FIG. 2: The predictions of the CQSM for A^ d (0) and A^ d (0) as functions of (the filled 
diamonds), together with the corresponding lattice predictions. Here, the open triangles corre- 
spond to the predictions of the LHPC group 18], while the open squares to those of the QCDSF 
collaboration ,21]. 



FigH shows the predictions of the CQSM for ^4^(0) and A$ d (0) as functions of m^, to- 
gether with the corresponding lattice predictions. As for AiQ d (0), the CQSM predictions and 
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the lattice QCD predictions are both independent of and consistent with the constraint 
of the quark number sum rule : 

A^ d (0) = N u + N d = 3, (127) 

with high numerical precision. Turning to A^ d (0), one finds a sizable difference between 
the predictions of the CQSM and of the lattice QCD. The lattice QCD predicts that 

A u + d {Q) = (x) u + (x) d ~ (0.6 ~ 0.7), (128) 

which means that only about (60 ~ 70)% of the total nucleon momentum is carried by the 
quark fields, while the rest is borne by the gluon fields. On the other hand, the CQSM 
predictions for the same quantity is 

A%+ d (Q) = (x) u + (x) d = 1, (129) 

which means that the quark fields saturates the total nucleon momentum. This may certainly 
be a limitation of an effective quark model, which contains no explicit gluon fields. Note, 
however, that the total quark momentum fraction A2Q d (0) is a scale dependent quantity. 
The lattice result corresponds to the energy scale of Q 2 = (2GeV) 2 while the CQSM 
prediction should be taken as that of the model energy scale around Q 2 = 0.30 GeV 2 ~ 
(560 Me V) 2 [40] . We shall later make more meaningful comparison by taking care of the 
scale dependencies of relevant observables. 

Next, in FigOl we show the isovector combination of the generalized form factors AiQ d (0) 
and A^ d (0). The meaning of the symbols are the same as in FigH As for A^ d (0), both 
the CQSM and the lattice simulation reproduce the quark number sum rule 

A^ d {0) = N u — N d = 1, (130) 

with good prediction. Turning to A%Q d (0), one observes that the prediction of the CQSM 
shows somewhat peculiar dependence on the pion mass. Starting from a fairly small value 
in the chiral limit (m n = 0), it first increases as m n increases, but as m n further increases 
it begins to decrease, thereby showing a tendency to match the lattice prediction in the 
heavy pion region. Very interestingly, letting put aside the absolute value, a similar 
dependence is also observed in the chiral extrapolation of the lattice prediction for the 
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FIG. 3: The predictions of the CQSM for (0) and A% d (0) as functions of m w , together with 
the corresponding lattice predictions ^^,[2^. The meaning of the symbols are the same as in 
FigEJ 



momentum fraction (x) u — (x) d shown in Fig. 25 of 17 1. Physically, the quantity d (0) 
has a meaning of the difference of the momentum fractions carried by the w-quark and the 
d-quark. The empirical value for it is A^ d (0) = (x) u ~ d = 0.154 ± 0.003 [rjj. One sees that 
the prediction of the CQSM in the chiral limit is not far from this empirical information, 
although more serious comparison must take account of the scale dependence of (x) u ~ d . 

Next, shown in Figgjare the CQSM predictions for B$ d (0) and B$ d (0). The former 
quantity is related to the isoscalar combination of the nucleon anomalous magnetic moment 
as B^ d (0) = k u + n d = 3(k p + n n ) = 3k (t=0) . (We recall that its empirical value is 
B^ d (0) ~ -0.36.) We find that this quantity is very sensitive to the variation of the 
pion mass. It appears that the CQSM prediction BiQ d (0) ~ —1.5 corresponding to chiral 
limit underestimates the observation significantly However, the difference is exaggerated 
too much in this comparison. In fact, if we carry out a comparison in the total isoscalar 
magnetic moment of the nucleon \G^~iq(0) = fi p + /i n = fi^ T=0 \ the CQSM in the chiral 
limit gives /J>cqsm — 0-5 i n comparison with the observed value /i^r - 1 ~ 0.88. To our 
knowledge, no theoretical predictions are given for this quantity by either of the LHPC or 
QCDSF collaborations. The right panel of Fig.4 shows the predictions for B%Q d (Q), which is 
sometimes called the isoscalar part of the nucleon anomalous gravitomagnetic moment, or 
alternatively the net quark contribution to the nucleon anomalous gravitomagnetic moment. 
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FIG. 4: The predictions of the CQSM for B^ d (0) and B%£ d (0) as functions of m w , together with 
the corresponding lattice predictions 19|.|22||. The meaning of the symbols are the same as in 
FigEJ 



As already pointed out, the prediction of the CQSM for this quantity is exactly zero : i.e. 



B£ d (0) = 0, 



;m) 



The explicit numerical calculation also confirms it. It should be recognized that the above 
result B^ d (0) = obtained in the CQSM is just a necessary consequence of the momentum 
sum rule and the total nucleon spin sum rule, both of which are saturated by the quark field 
only in the CQSM as 

A^ d (0) = (x) u+d = 1, (132) 



and 



~ K d (o) + B>£ d (o)} = (jy+ d = i. 



(133) 



In real QCD, the gluon also contributes to these sum rules, thereby leading to more general 
identities : 

^o +d (0) + ^ (0) = l, (134) 
[A^ d (0) + B^ d (0)] + [A%(0) + Sf (0)] = 1, (135) 

which constrains that only the sum of B%Q d (0) and -Bfo(O) ^ s f° rce d to vanish as 

B£ d (0) + Bf (0) = 0. (136) 
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(While we neglect here the contributions of other quarks than the u- and <i-quarks, it loses 
no generality in our discussion below. In fact, to include them, we have only to replace the 
combination u + d by u + d + s + ■ ■ ■.) The above nontrivial identity claims that the net 
contributions of quark and gluon fields to the anomalous gravitomagnetic moment of the 
nucleon must be zero. An interesting question is whether the quark and gluon contribution 
to the anomalous gravitomagnetic moment vanishes separately or they are both large with 
opposite sign. A perturbative analysis based on a very simple toy model indicates the latter 
possibility |32| . On the other hand, a nonperturbative analysis within the framework of the 
lattice QCD indicates that the net quark contribution to the anomalous gravitomagnetic 
moment is small or nearly zero, i^o^O) = Q],!^]- (To be more precise, we sees that 
the prediction of the LHPC collaboration for B%Q d (Q) is slightly negative Q], while that of 
the QCDSF group is slightly positive [22].) This strongly indicates a surprising possibility 
that the quark and gluon contribution to the anomalous gravitomagnetic moment of the 
nucleon may separately vanish. Worthy of special mention here is an interesting argument 
given by Teryaev some years ago, claiming that the vanishing net quark contributions to 
the anomalous gravitomagnetic moment of the nucleon, violated in perturbation theory, is 
expected to be restored in full nonperturbative QCD due to the confinement j^, 48^ 49^ . 
Very interestingly, once it actually happens, it leads to a surprisingly simple result, i.e. the 
proportionality of the quark momentum and angular momentum fraction 

J u+d = i (x) u+d , (137) 

as advocated by Teryaev j^], 48^!3- A far reaching physical consequence resulting from 
this observation was extensively discussed in our recent report [46]. (See also the discussion 
at the end of this section.) 

Next, we show in FigJSJthe predictions for the isovector case, i.e. B^Q d (0) and B%Q d (0). 
We recall first that the quantity 5"^" d (0) represents the isovector combination of the nucleon 
anomalous magnetic moment k^ t=1 ^ = k p — n n ( = k u — n d = k^ i=1 ^), the empirical value of 
which is known to be k^ t=1 ^ = 3.706. One find that this quantity is extremely sensitive to 
the variation of the pion mass especially near m n = 0. This is only natural if one remembers 
the important role of the pion cloud in the isovector magnetic moment of the nucleon. (One 
may notice that the prediction of the CQSM for k^ t=1 ^ underestimates a little its empirical 
value even in the chiral limit. We recall, however, that, within the framework of the CQSM, 
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FIG. 5: The predictions of the CQSM for B^ d (0) and B^ d (ti) as functions of m w , together with 
the corresponding lattice predictions 19] ,'22]. The meaning of the symbols are the same as in 
FigEJ 

there is an important 1/N C correction or the 1st order rotational correction to some kind 
of isovector quantities like the isovector magnetic moment of the nucleon in question or the 
axial-vector coupling constant of the nucleon This next-to-leading correction in 

1/N C should also be taken into account in more advanced investigations.) Shown in the 
right panel of Fig|5] is the theoretical predictions for B%Q d (0), the half of which can be 
interpreted as the difference of the total angular momentum carried by the w-quark and the 
(i-quark fields according to Ji's angular momentum sum rule j^j]. The CQSM predicts fairly 
small value for this quantity, in contrast to the lattice predictions of sizable magnitude. 
It seems that the pion mass dependence rescues this discrepancy only partially. Here we 
argue that, the reason why the CQSM (in the chiral limit) gives rather small prediction for 
this quantity is intimately connected with the characteristic x dependence of the quantity 
E^m~ 1 \%, 0, 0), the forward limit of the isovector unpolarized spin-flip GPD of the nucleon. 
To show it, we first recall that, within the theoretical frame work of the CQSM, BiQ d (Q) as 
well as B%Q d (0) are calculated as difference of G^-^(O) and G^q\o) and of G^f^(0) an d 
(j£20 (0)> respectively, as 

B£*(P) = G^(0) - G£ 1} (0), (138) 
B u 20 d (0) = 4/3(0) - Ggo 1} (0). (139) 



Although the quantities of the rhs can be calculated directly without recourse to any dis- 
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tribution functions, they can also be evaluated as x-weighted integrals of the corresponding 
GPDs as 

(140) 
(141) 

J 1 f u ~ d {x) dx = N u -N d = 1, (142) 
J xf u - d {x) dx = (x) u - (x) d . (143) 

The distribution function Ej^~ 1) (x, 0, 0) has already been calculated within the CQSM in 
our recent paper (2^. As shown there, the Dirac sea contribution to this quantity has a 
sizably large peak around x — 0. Since this significant peak due to the deformed Dirac- 
sea quarks is approximately symmetric with respect to the reflection x — > —x, it hardly 
contributes to the second moment G^^(Q), whereas it gives a sizable contribution to the 
first moment Gjjfio (0). The predicted significant peak of E^f (x, 0, 0) around x = can 
physically be interpreted as the effects of pion cloud. It can be convinced in several ways. 
First, we investigate how this behavior of (x, 0, 0) changes as the pion mass is varied. 
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FIG. 6: The prediction of the CQSM for E%f d (x, 0,0) = ^£ =1) (x, 0, 0) obtained with M = 
400 Me V and m n = 0. 

Shown in FiglHl and in Fig|7| are the CQSM predictions for E^~ X \x, 0, 0) with several 
values of m n . i.e. m n = 0,200, and 400 MeV. One clearly sees that the height of the peak 



G ( m%\0) = £ E i I !r 1 \x,0,0)dx, 

GgS(0) = f i xE{ I r 1 \x,0,0)dx, 

Gg } (0) = j' i H i J =1 \x J 0,0)dx = 

&ijo\ Q ) = J l xH { E =1 \x,0,0)dx = 
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FIG. 7: The m n dependence of E%f d (x, 0, 0). 

around x — 0, due to the deformed Dirac-sea quarks, decreases rapidly as increases. This 
supports our interpretation of this peak as the effects of pion clouds. On the other hand, one 
also observes that the magnitude of the valence quark contribution, peaked around x ~ 1/3, 
gradually increases as becomes large. This behavior of E^ [ ~ 1 \x, 0, 0) turns out to cause 
a somewhat unexpected dependence of G^~^(0) and G^^(0). As a function of m n , the 
Dirac sea contribution to G^q^O) decreases fast, whereas the valence quark contribution 
to it increases slowly, so that the total G^io^O) becomes a decreasing function of m n . On 
the other hand, owing to the approximate odd-function nature of the Dirac sea contribution 
to xE <y l [ ~ l \x ) 0, 0) with respect to x, it hardly contributes to G^f^iO) independent of the 
pion mass, while the valence quark contribution to xE^ 4 ~ 1 \x, 0, 0) is an increasing function 
of m^, thereby leading to the result that the net G^^(0) is a increasing function of m n . 



We can give still another support to the above-mentioned interpretation of the contri- 
bution of the Dirac-sea quarks. To see it, we first recall that the theoretical unpolarized 
distribution function f u ~ d (x) appearing in the decomposition 

E%f d {x, 0, 0) = f u ~ d (x) + E u - d (x, 0, 0), (144) 

also has a sizable peak around x = due to the deformed Dirac-sea quarks. As shown in 
Fig |H] and in Fig El this peak is again a rapidly decreasing function of m n , supporting our 
interpretation of it as the effects of pion clouds. Here, we can say more. We point out that 
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X 

FIG. 8: The prediction of the CQSM for f u - d (x,0,0) obtained with M = 400 MeV and = 0. 




FIG. 9: The m w dependence of f u - d (x,0,0). 

;his small- x behavior of f u ~ d (x) is just what is required by the famous NMC measurement 
53^ . To confirm it, first remember that the distribution f u ~ d (x) in the negative x region 
should actually be interpreted as the distribution of antiquarks. To be explicit, it holds that 

u{x) - d(x) = -f u ~ d (-x) (x > 0). (145) 



The large and positive value of f u (x) in the negative x region close to x = means that 

c ove 
53]. 



u(x) — d(x) is negative, i.e. the dominance of the d-quark over the w-quark inside the proton, 
which has been established by the NMC measurement 
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FIG. 10: The prediction of the CQSM for d(x) - u(x) evolved to Q 2 = 2.3 GeV 2 and Q 2 = 54 GeV 2 
in comparison with the Hermes and NuTeV data at the corresponding energy scales |57l |. 5Jj|. 



Shown in Fig^Jare the predictions of the CQSM for d(x)—u(x) evolved to the high energy 
scales corresponding to the experimental observation 54] . (The theoretical predictions here 
were obtained with M = 400 MeV and = 100 Me V.) The model reproduces well the 
observed behavior of d(x) — u(x), although the magnitude of the flavor asymmetry in smaller 
x region seems to be slightly overestimated. It is a widely accepted fact that this flavor 
asymmetry of the sea quark distribution in the proton can physically be understood as the 
effects of pion cloud at least qualitatively 3|-44|. This then supports our interpretation 
of the effects of the deformed Dirac-sea quarks in E ( ^~ 1 \x, 0, 0) and f u ~ d (x) as the effects 
of pion clouds. 

We show in Table [Q the model predictions for G^{0), G^aJ^O), and B^ d (0) = 



Gmw (°) 



G^ 2 q^(0) as functions of m^. One sees that the value of G^ 2 ^(0) with m n = 



is an order of 0.3 ~ 0.4. As already pointed out, it gradually increases as becomes large. 
This is also the case with G { ew\Q)- As 

a consequence, the isovector combination of the 
nucleon anomalous gravitomagnetic moment i?2o~ d (0), which is obtained as a difference of 
the above two quantities, is also an increasing function of m^, thereby having a tendency 
to come closer to the lattice prediction given in the heavy pion region. Still, the CQSM 
predictions B^ d (0) ~ 0.3 around = 500 MeV is a factor of two smaller than the corre- 
sponding lattice prediction i?2o~ d (0) — 0.6. Now we summarize the reason why the CQSM 
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TABLE I: The dependencies of G$£j(0), G%^\o), and B% d {0) in the CQSM with M = 
400 MeV. 



(MeV) 


u M,20 V / 


G^on^O) 
E,20 \ ' 


20 v / 


n 
u 


u.OUl 


998 


U. lOO 


100 


0.392 


0.276 


0.116 


200 


0.452 


0.327 


0.125 


300 


0.519 


0.350 


0.169 


400 


0.579 


0.354 


0.225 


500 


0.640 


0.347 


0.293 


600 


0.716 


0.328 


0.388 



gives fairly small prediction for B%q (0). It is due to two types of cancellations. The first 
is the cancellation of the potentially large contribution of Dirac-sea quarks arising from the 
approximately antisymmetric behavior of the Dirac sea contribution to xE^f 1 \x 1 0, 0) as 
well as xf u ~ d (x). The second is the cancellation between the total gravitomagnetic moment 
Gm~2q(0) and its canonical part G^q\o). We are not sure whether the lattice simulation 
carried out in the heavy pion region with neglect of the so-called disconnected diagrams can 
efficiently take account of such effects of chiral dynamics as discussed above. 

So far, we have investigated the pion mass dependence of the forward limit of the general- 
ized form factors of the nucleon. Here, we investigate the momentum-transfer dependencies 
of some form factors of the nucleon. For the reason explained before, all the physical pre- 
dictions given hereafter will be obtained with use of the mass parameter M = 400 MeV and 
= 100 MeV. We show in FiglTTl FiglT^l Fig fEfl Fig^^the predicted momentum-transfer 
dependencies of the generalized Sachs form factors, G^Yo 20 (0> ^bIo ^£20 (0> 

^M,10 (^)' ^M,20 (X)> ^M,10 ft)' an d ^M,20 (^)- 

To get some feeling about the momentum-transfer dependencies of the predicted form 
factors, we shall compare them with the existing empirical data. At present, only the lowest 
moment of the GPDs, i.e. the standard electromagnetic form factors of the nucleon are 
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FIG. 11: The predictions of the CQSM for the isoscalar generalized electric form factors ±*q (t) 
and 4 J = 0) (i). 
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FIG. 12: The predictions of the CQSM for the isovector generalized electric form factors G% iq (t) 
and G% 2c! (*)• 



experimentally known. Shown in FigEHare the predictions of the CQSM for the Dirac form 
factors of the proton and the neutron, in comparison with the empirical data 159] together 

nn 

with the corresponding predictions at the LHPC lattice simulation [18J, 19J. One observes 
that the t-dependence of the CQSM prediction for Ff (t) is a little too stronger than the 
empirical one, while the t-dependence of the lattice predictions is too much weaker than 
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FIG. 13: The predictions of the CQSM for the isoscalar generalized magnetic form factors G^ $ (t) 
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FIG. 14: The predictions of the CQSM for the isovector generalized magnetic form factors iq (t) 
and G { ^{t). 

the empirical one. A little too fast falloff of the CQSM predictions means that it slightly 
overestimates the electromagnetic size of the proton. On the contrary, the electromagnetic 
proton size predicted by the LHPC simulation is too small as compared with the empirically 
known size. As is well known, the Dirac form factor of the neutron is not well determined 
experimentally. Both of the CQSM prediction and the lattice QCD predictions are very 
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small in magnitude in qualitatively consistent with the empirical information, although the 
former is slightly positive, while the latter is slightly negative. 
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FIG. 15: The prediction of the CQSM for the Dirac form factors of the proton and the neutron in 
comparison with the empirical data 59] together with the LHPC lattice predictions. The LHPC 
lattice predictions are based on their dipole fits 

FiglTT)! shows the predictions of the CQSM for the Pauli form factors of the proton and the 
neutron, in comparison with the empirical data together with the corresponding predictions 
of the LHPC lattice simulation. (Here, both of the CQSM predictions and the lattice QCD 
predictions are normalized to the observed anomalous magnetic moments of the proton and 
the neutron at the zero momentum transfer.) The solid curves represent the predictions of 
the CQSM, whereas the dashed curves show the corresponding lattice predictions. One can 
see that the predictions of the CQSM reproduce the empirical Pauli form factors fairly well. 
On the other hand, the lattice QCD predicts too slow falloff of the Pauli form factors, which 
means that the magnetic sizes of the nucleon are largely underestimated by the lattice QCD 
simulations. The underestimate of the nucleon electromagnetic sizes seems to be a general 
tendency of the lattice QCD simulations by the LHPC and QCDSF collaborations. It is not 
clear yet whether the origin of discrepancy can be traced back to the fact that these lattice 
simulations were carried out in the region with unrealistically heavy pion mass. 

Concerning the genuine generalized form factors with n > 2, we have no experimental 
information yet. Among them, of our particular interest is the generalized form factor 
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FIG. 16: The predictions of the CQSM for the Pauli form factors of the proton and the neutron in 
comparison with the empirical data 
dipole fits 0,0- 



59( together with the LHPC lattice predictions based on the 



#20" (t) appearing in Ji's angular momentum sum rule. We show in Figll7lthe prediction 
of the CQSM for B^ d (t) obtained with M = 400 MeV and m w = 100 MeV, in comparison 
with the corresponding predictions of the LHPC group. One sees that, for arbitrary value of 
t, B2Q d {t) does not vanish in the CQSM, which is an indication of the fact that the shapes of 
the quark momentum and the total angular momentum distributions are not completely the 
same. Still, the magnitude of B2Q d (t) turns out to be very small, which seems qualitatively 
compatible with the prediction of the LHPC group, although one should not forget about 
large uncertainties in the lattice simulation at the present level. 

We are now ready to discuss what we can say about the spin contents of the nucleon 
from our investigation on the generalized form factors of the nucleon. The quantities of our 
interest are all obtained from the forward limit of the generalized form factors, which are 
denned as the second moments of the relevant GPDs : 



(x) u+d 
2J u+d 



G 
G 



(1=0) 
E,20 

(1=0) 
M,20 



(0), 
(0), 



2 J u ~ 



(1=1) 
M,20 



G 



(0), 
(0). 



(146) 
(147) 



Summarized below are the predictions of the CQSM model for these quantities obtained 



35 



0.5 



0.0 



-0.5 



< 


B 2 r (t) 

^ i — ^- 




i.i. 



0.0 



1.0 



2.0 



3.0 



-t [GeV 2 ; 



FIG. 17: The predictions of the CQSM for the generalized form factor B2Q d (t) of the nucleon, in 
comparison with the corresponding prediction of the lattice QCD by LHPC collaboration [ig| ] . 



with M = 400 MeV and m„ 



100 Me V : 

u+d 



(x) 
2J u+d 



1.00, 
1.00, 



(x) u ~ d = 0.276, 

2 JU-d = Q_ 406 _ 



(148) 
(149) 



As pointed out before, these second moments of GPDs are generally scale dependent. Our 
viewpoint is that the predictions of the CQSM correspond to those at the low energy scale 
where the validity of the model is ensured. (This energy may typically be characterized by 
the Pauli-Villars mass Ai ~ 600 MeV.) We shall take account of the scale dependence of the 
above quantities by solving the QCD evolution equation at the next-to-leading order (NLO) 

nn 

with the predictions of the CQSM as the initial conditions 60j-[62]. For simplicity, let us 
assume that, at this low energy scale, there is no contribution of gluon fields as well as of 
the strange quarks, which dictates that 

(x) s = 0.0, 2 J s = 0.0, (x) 9 = 0.0, 2J 9 = 0.0. (150) 

The starting energy of the evolution is taken to be Q\ ni = 0.30 GeV 2 ~ (550 MeV) 2 , because 
it is favored from the previous successful application of the model to high energy deep- 
inelastic-scattering observables Q-Q- Taking N f = 3 and A QCD = 0.248 GeV, we find 
that, at Q 2 = 4 GeV 2 , 

(x) u+d+s = 0.676, (x) u ~ d = 0.171, (x) 9 = 0.324, (151) 
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2JU+d+s = q 6?6) 2JU - d = q 2jg = q 



(152) 



One may notice that the values of (x) u+d+s and 2 J u+d+s at Q 2 = 4 GeV 2 precisely coincide. 
Actually, the equality of these two quantities holds at any energy scale. The reason is because 
these two quantities obey exactly the same evolution equation and because they are equal 



atter is also 

2210 In 
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at the initial energy scale according to the CQSM jllj. (We emphasize that the 
the case for the LHPC and QCDSF lattice predictions at least approximately 
table HH we compare these predictions of the CQSM with those of the lattice QCD and also 
with the empirical values for (x) u+d+s and (x) u ~ d obtained from the phenomenological PDF 
fits 64]. As one sees, the total momentum fraction (x) u+d+s carried by the quarks decreases 
rapidly as Q 2 increases. The evolved value (x) u+d ~ 0.68 at Q 2 = 4 GeV 2 is not extremely 
far from the lattice QCD prediction at the same normalization scale, although it is a little 
larger than the empirical value (x)emptricai — 0.57. This difference may be an indication of 
the fact that, even at the low energy scale around Q 2 ~ 0.30 GeV 2 , the gluons may carry 
some portion of the nucleon momentum. As far as the difference of the momentum fractions 
carried by the u-quark and the (i-quark, the CQSM well reproduces the empirical value 
i x ) empirical — 016 at Q 2 = 4 GeV 2 , whereas the lattice QCD overestimates it a little. 

TABLE II: The predictions of the CQSM for 2 J u+d , and 2 J u ~ d in comparison 

with the predictions of the lattice QCD simulations Q| , Q] , [3] , 01 as well as with the empirical 
information jf^j]. 





CQSM (model scale) 


CQSM (Q 2 = 4 GeV 2 ) 


LHPC 


QCDSF 


empirical 


(x) u+d+s 


1.000 


0.676 


0.61 


0.59 


0.57 


(x) u ~ d 


0.276 


0.171 


0.269 


0.24 


0.157 


2 ju+d+s 


1.000 


0.676 


0.58 


0.66 




2J u-d 


0.406 


0.257 


0.93 


0.82 





Next we turn to the discussion of the total angular momenta J u and J d carried by the 
quark fields, on which we do not have any empirical information yet. One can see that, as far 
as the total angular momentum fraction J u+d+s carried by the quark fields, is concerned, the 
prediction of the CQSM is qualitatively consistent with that of the lattice QCD. However, 
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a big discrepancy is observed for the difference J u ~ d of the angular momentum carried by 
the w-quark and the (i-quark. The cause of this discrepancy can be traced back to that 
of the isovector gravitomagnetic moment of the nucleon G^^(O), which we have already 
discussed. Fairly small prediction of the CQSM for J u ~ d also appears to be incompatible 
with the semi-theoretical (or semi-phenomenological) estimate carried out in [6| with partial 
use of the predictions of the CQSM. As we shall discuss below, however, their estimate for 
J u and J d shown in Table 4 and Table 5 of [6] should be taken with care. In fact, it was 
obtained based on the valence-like approximation, i.e. by neglecting the sizable Dirac sea 
contributions to E u (x,0,0) and E d (x,0,0). To be more concrete, they start with a simple 
guess for these distributions functions as 

E u (x,0,0) = ~K u f u ™>(x), (153) 

E d {x,0,0) = K d f d ™ l {x), (154) 

with 

K u = 2k p + k u = 1.673, (155) 

K d = K p + 2k u = -2.033. (156) 

This parameterization trivially satisfies the 1st moment sum rule 

J E q (x,0,0)dx = n q . (157) 

On the other hand, by using the 2nd moment sum rule or Ji's angular momentum sum rule, 
they obtain 

J u = ~ [(x} u + k u (x) u ™'} } (158) 



J" 



1 r 



(x) d + n d (x) aval , (159) 



where (x) q is the momentum fraction carried by the quark of flavor q, while (x) qval is the 
corresponding contribution of the valence quark in the sense of parton model. Using the 



MRST98 parameterization for the unpolarized PDFs [631] . they could thus obtain at Q ~ 
lGeV 2 

(x) Uval ~ 0.34, (x) u ~ 0.40, 2 J u ~ 0.69, (160) 

and 

(x) dval ~ 0.14, (x) d ~0.22, 2J d ~-0.07, (161) 
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which appears to be qualitatively consistent with the lattice predictions. (Here, we have 
discarded the small contribution of the s-quark, for simplicity.) However, after this simple 
estimate for J u and J d , they next try to take account of the sizable Dirac sea contribution to 
the distributions E u (x,0,0) and E d (x 1 0, 0). As already shown in our exact model calcula- 
tion, and as shown in Fig.8 of which is obtained based on the derivative-expansions-type 
approximation within the CQSM, the Dirac sea contribution to E u {x, 0, 0) has a narrow 
and positive peak around x ~ 0. To simulate this narrowly peaked behavior of the Dirac-sea 
contribution, they propose to parameterize it by a 5 function in x. The new and improved 
parameterization for E u (x, 0, 0) and E d (x, 0, 0) are then given as 

E u {x,0,0) = A u f u ^(x) + B u S(x), (162) 
E d {x,0,0) = A d f d ^{x)+B d 5{x). (163) 

From the 1st and 2nd sum rules for E u (x, 0, 0) and E d (x, 0, 0), they obtain 



9 T u — ( T ) 1 



A* = 2 ■*-<*>', (165) 

(x} dval 



B u 



1 „ 2 J u - {x) u 

2 (x) Uml 

2 J d - (x ) d 
(x) dva 



(166) 



B d = K d - x ,7 ■ (167) 



As pointed out in |6(, the total angular momentum carried by u- and <i-quarks, J u and J d , 
now enter as fit parameters in the parameterization of E u (x,0,0) and E d (x,0,0). If so, 
there is no compelling reason to believe that they are close to the estimate given in ()158j) 
and ()159j) . obtained within the valence-type parameterization for E u (x, 0, 0) and E d (x, 0, 0). 
In fact, if one puts the estimate given in (|16U|) and (|161|) into the above relations (|165j) ~ 
(|167j) . one obtains 

A u ~ 0.853, A d ~ -2.071, B u ~ -0.033, B d ~ 0.038, (168) 

which dictates that the coefficients of the 5 functions are small or nearly zero. This is only 
natural, since the used value of J u and J d are just estimated based on the valence- type 
parameterization for E u (x,0,0) and E (x, 0,0). Conversely speaking, the values of J u and 
J d quoted in Table 4 and table 5 of p] need a revision, because they are incompatible with 
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the existence of sharp peak of E u d (x,0, 0) with sizable positive magnitude around x ~ 
observed in Fig.8 of the same paper. The inseparable relation between the magnitude of J u — 
J d and the sharp of E u ~ d (x, 0, 0) can be made more transparent by slightly modifying their 
schematic analysis. Instead of E u ~ d (x, 0, 0), we chose here to parameterize E^f d (x, 0, 0) = 
H u ~ d (x, 0, 0) + E u - d (x, 0, 0) as 

E u M - d (x,0,0) = ci [f u ™ l (x)-f d ™ l (x)] + c 2 5(x). (169) 

Here the 5 function term is thought to simulate the sizable sharp peak of E~^f d (x, 0, 0) 
predicted by the CQSM. Actually, it need not be a 5 function. It can be an any function g 
of far as it satisfies the following two conditions : 

• g(x) is an even function of x at least approximately, 

• the integral of g(x) over x gives a (positive) number c 2 . 

Assuming that these conditions are satisfied (as is the case for the predictions of the 
CQSM), the 1st and the 2nd moment sum rule of E'^j~ d (x, 0, 0) leads to the identities : 

1 + K u -n d = J 1 E%f d (x, 0, 0) dx = a + c 2 , (170) 

2 (J u - J d ) = J 1 El[ d {x, 0, 0) dx = ci ( (x) u ™ 1 - (x) d ™ 1 ) . (171) 

Combining these relations, we have 

2 (J u - J d ) = i ( - (x)^^ ) . (172) 

Here we have introduced a parameter r = c 2 /ci. This relation shows that, except for one 
parameter r, the difference of the total angular momenta, carried by the u- and <i-quarks, 
ju _ jd^ j g gj ven by K d^ ( x }uvai anc i (^x) dva \ which are all observables. What is the physical 
meaning of the parameter r, then ? To understand it, we first recall that the quantity 
E^ d {x ) 0, 0) represents the distribution of the isovector magnetic moment of the nucleon 
in Feynman x-space not in ordinary coordinate space. The contribution of its valence-like 
distribution to the magnetic moment gives c±, whereas that of its sea-like distribution gives 
c 2 . The CQSM indicates that they are of approximately equal magnitude, i.e. c\ ~ c 2 , which 
means that r ~ 1. On the other hand, roughly speaking, the lattice simulations carried out 
in the heavy pion region with neglect of the disconnected diagrams correspond to r ~ 0. 
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TABLE III: The value of 2 (J u - J d ) functions of the ratio of the parameters C2 and c\ defined 
in (HI]) . 



r = c 2 /ci 


2 (J u - J d ) 





0.941 


0.5 


0.627 


1 


0.471 



Table ILTT1 show the values of 2(J U — J d ) obtained from (|172|) for several typical values of 
the ratio r. The quantity (x) Uval — (x) dval are actually scale dependent. For simplicity, here 
we have used the empirical value (x) Uval — (x) dval ~ 0.2 corresponding to Q 2 = 1 GeV 2 quoted 
in 0| . One sees that the two cases, i.e. r = case and r = 1 case, lead to a factor of two 
difference for J u — J d . What is indicated by this observation is an inseparable connection 
between the angular momentum carried by the quark fields in the nucleon and the magnetic 
moment of the nucleon, or more precisely the distribution of magnetic moments in Feynman 
x-space. As a matter of course, the relation between the quark angular momentum and the 
nucleon magnetic moment could be anticipated from more general ground of Ji's angular 
momentum sum rule. However, an advantage of our explicit model analysis is that we can 
get more deep and concrete insight into the possible behavior of the relevant distribution 
E u ~ d (x, 0, 0), on which we have no experimental information yet. 

Returning to the isoscalar (flavor-singlet) combination or the net quark contribution to 
the total angular momentum of the nucleon, J u+d+s ; we have pointed out that the predic- 
tion of the CQSM for it is not so far from that of the lattice QCD. However, both of LHPC 
and QCDSF collaborations also estimated the net orbital angular momentum carried by the 
quark fields, thereby being led to the conclusion that the total orbital angular momentum 
of the quarks is very small or consistent with zero. As pointed out in our recent paper, this 
conclusion contradicts not only the prediction of the CQSM but also the famous EMC ob- 
servation. Although the possible reason of this discrepancy was already pointed out in that 
paper [46]. here we discuss it in more detail, especially by taking care of the scale dependen- 
cies of the relevant observables. The quark orbital angular momentum can be obtained by 
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subtracting the intrinsic quark spin term from the total quark angular momentum J u+d as 

L u+d = ju+d _ 1 AS «+d ? ^73) 

where J u+d is given by 

J u+d = \ ( (x) u+d + B^ d (0)) , (174) 

Using the results of the dipole fits for the generalized form factors, the LHPC collaboration 
obtain 

( x ) u + d = A^ d (0) = 0.666 ±0.009, (175) 
££ (0) = 0.29 ±0.04, B d (0) = -0.38 ±0.02, (176) 
AS" = 0.860 ±0.069, AS d = - 0.171 ± 0.043, (177) 

which in turn gives 

L u = -0.088 ±0.019, L d = 0.036 ± 0.013, (178) 

or 

L u+d = _ 0.052 ±0.019. (179) 
On the other hand, the QCDSF collaboration obtain 

(x) u = A% (0) = 0.400 ±0.022, (x) d = A d (0) = 0.147 ±0.011, (180) 

B% (0) = 0.334 ±0.113, B d (0) = - 0.232 ± 0.077, (181) 

AS" = 0.84 ±0.02, AS d = -0.24 ±0.02, (182) 

which gives 

L u+d = o.03 ±0.07. (183) 

As one sees, a common conclusion of the two groups is that the total orbital angular mo- 
mentum of quarks is very small or consistent with zero. 

Since these lattice predictions corresponds to the energy scale of Q 2 = 4 GeV 2 in the MS 
scheme, we try to evolve the corresponding predictions of the CQSM to the same energy. 
To know the scale dependence of the total quark orbital angular momentum L Q with Q 
denoting the sum of all quark flavors, we need to know the scale dependence of the total 
quark angular momentum J Q and that of the total quark longitudinal polarization AS Q . 
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We recall again the fact that the angular momentum fractions carried by the quark and the 
gluon field, i.e. J® and J 9 , obey exactly the same evolution equation as the total momentum 
fractions carried by the quark and gluon fields (x)® and (x) 9 The evolution equation 
of A£^, which is coupled with the evolution of the gluon polarization, is also well known 
69^ . As initial conditions of the evolution, we use the predictions of the CQSM (the flavor 
SU(2) version) : 



2J 



u+d 



1.0, 



;i84) 



supplemented with the assumption 



2 J s = 0.0, A£ s = 0.0, Ag = 0.0, 



(185) 



421. 



at Qf ni = 0.30 GeV 2 . (There also exists the flavor SU(3) version of the CQSM |41| 
It predicts that A£ s is a negative quantity of the order of (5 ~ 10) %. However, the 
flavor singlet combination or the net quark contribution to the total quark longitudinal 
polarization AS^ takes almost the same value in both version of the CQSM, so that the 
following discussion will receive no modification.) 




2.0 3.0 
Q 2 [GeV 2 ] 




FIG. 18: The left panel shows the scale dependencies of the quark and gluon total angular momenta 
J® and J 9 obtained by solving the evolution equation at the NLO with the predictions of the CQSM 
as initial conditions, while the right panel shows the scale dependence of A£^ and Ag. 



The left panel of FiglTHl shows the scale dependence of J® obtained by solving the evolu- 
tion equation at the NLO in the MS scheme. (Here we set N c = 3 and Aq C d = 0.248 GeV.) 

43 



One observers that, especially at low energy scales, J Q is a rapidly decreasing function, 
while J 9 is a rapidly increasing function of Q 2 . The right panel of FigfTBI shows the scale 
dependence of AE^ and Ag obtained by solving the NLO evolution equation in the same 
renormalization scheme As one sees, Ag is a rapidly- increasing function of Q 2 . On the 
other hand, AE u+d+s has a fairly weak scale dependence. Its scale dependence is restricted 
to the very low energy region below Q 2 < 0.6 GeV 2 and beyond that scale it changes very 
slowly. (We recall that, at the leading order (LO), AS^ in the MS is exactly scale inde- 
pendent. One may also remember the fact that AE^ in the chiral-invariant renormalization 
scheme is scale independent by definition .) Now, combining the results for the scale 

dependencies of J Q , J 9 , AE Q , and Ag, one can predict the scale dependencies of L Q and L 9 
at the NLO from 

L Q = J Q -^AZ Q , (186) 
L 9 = J 9 -Ag. (187) 

FiglT9lshows the scale dependencies of quark and gluon orbital angular momenta obtained in 
this manner. One sees that the total quark orbital angular momentum is a rapidly decreasing 
functions of the energy scale, especially at the low energy scales. Since the longitudinal quark 
polarization is only weakly scale dependent, this feature comes from the scale dependence 
of the total quark angular momentum, which has the same scale dependence as the total 
quark momentum fraction. One also sees that the gluon orbital angular momentum is a 
decreasing function of the energy scale. 

For the sake of comparison with the lattice QCD predictions corresponding the the energy 
scale of Q 2 = 4 GeV 2 , we summarize in Table ITVl the predictions of the CQSM for the 
nucleon spin contents at the same energy scale. One confirms that the total quark orbital 
angular momentum is a rapidly decreasing function of the energy scale and its value at 
Q 2 = 4 GeV 2 is nearly half of that at the low energy model scale around Q 2 ni = 0.30 GeV 2 . 
Nonetheless, it still bears a sizable amount of the total nucleon spin even at the scale 
Q 2 = 4 GeV 2 , in contrast to the lattice predictions. We point out that, after taking account 
of the scale dependence, the predictions for J® are not so different between the CQSM and 
the lattice QCD. What is remarkably different is the predictions of the two theories for the 
net longitudinal quark polarization or the contribution of intrinsic quark spin. It is clear 
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FIG. 19: The scale dependencies of the quark and gluon orbital angular momenta and L 9 at 
the NLO obtained with combined use of the predictions of the CQSM and the QCD evolution 
equations at the NLO. 

that the lattice QCD simulations by the LHPC and the QCDSF collaborations for AS^ 
considerably overestimate the empirically known value of AE^, which is known to be quite 
small as 

Empirical = (0-2 ~ 0.35), (188) 

while the prediction of the CQSM is qualitatively consistent with this empirical information. 
A plausible reason why the lattice simulations by the LHPC and the QCDSF collaboration 
predicts fairly large AS Q around 0.6 was pointed out in 2^|. In that paper, we investigated 
the pion mass dependence of AE^ within the framework of the CQSM and found that 
it is very sensitive to the variation of m^, especially in the region close to the chiral limit 
nix = 0. (The magnitude of AS^ decreases rapidly as m w approaches 0.) This indicates that 
the lattice estimates carried out in the heavy pion region around m n = (700 ~ 900) MeV 
may not give reliable prediction for the particular observable AS^. As a consequence, their 
conclusion that the orbital angular momentum carried by the quark fields in the nucleon 
is negligible, must also be taken with care. It may be justified in the heavy pion world, 
but whether it is also the case in our chiral world is a different question, which should be 
answered by the lattice QCD studies in the future. 

Also noteworthy is the fact that the large values of AT, u+d obtained by the LHPC and 
QCDSF lattice collaborations seem to contradict the results of earlier lattice QCD studies 



45 



74|-[3|, which predict fairly small AS^ around (0.2 ~ 0.3). Among them, Mathur et. al. 

74j also estimated the total quark angular momentum J® from the quark energy-momentum 
tensor form factors on the lattice with the quenched approximation, and deduced that the 
quark orbital angular momentum carries about 34 % of the total proton spin, which is 
compatible or even dominant over the contribution of intrinsic quark spin around 26 % 
obtained by their simulation. 

TABLE IV: The predictions of the CQSM for the spin contents of the nucleon in comparison with 
the corresponding predictions of the LHPC and QCDSF lattice QCD simulations 
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CQSM (model scale) 


CQSM (Q 2 = 4GeV 2 ) 


LHPC 


QCDSF 


2 ju+d+s 


1.000 


0.676 


0.56 


0.66 




0.350 


0.318 


0.69 


0.60 


2 jja+d+s 


0.650 


0.358 


-0.11 


0.06 



So far, we have given several interesting theoretical predictions on the basis of an effective 
theory of QCD, i.e. the CQSM. Although we believe that the reliability of these predictions is 
guaranteed by the phenomenological success of the CQSM achieved in the nucleon structure 
function physics, it would be nicer if one can extract some predictions, which do not depend 
on a specific model of the nucleon. It is in fact possible, if one accepts the following two 
theoretical postulates. They are 

• Ji's angular momentum sum rule : J Q = \ [(x) Q + -£>2o(0)] , 

• absence of the net quark contribution to the anomalous gravitomagnetic moment of 
the nucleon : B$ (0) = . 

It is reasonable to accept the first postulate. Otherwise, we would lose a only clue to 
experimentally access the quark angular momentum in the nucleon. What is crucial in the 
following argument is therefore the second postulate. As already mentioned, the identity 
£?£q(0) = 0, that holds within the CQSM, just follows from the total momentum and the 
total angular momentum sum rules, both of which are saturated by the quark fields alone 
in this effective quark theory. It can therefore be an artifact of the model. However, an 
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interesting observation is that the smallness of -820(0) is also the predictions of the LHPC 
and the QCDSF lattice QCD simulations, which take account of full quark-gluon dynamics. 
Naturally, one should not forget about large uncertainties in the lattice simulations at the 
present stage. One should also worry about the m^-dependence of 5^,(0), although we 
conjecture from our analyses in the CQSM a weak m^-dependence of this quantity. 

Here, we shall proceed by assuming that -820(0) vanishes exactly or at least very small. As 
already pointed out, the identity -B^O) leads to remarkable relations, i.e. the proportionality 
of the total momentum and total angular momentum carried by the quark fields and also 
by the gluon fields as 

J Q = \(x) Q , J 9 = \(xY. (189) 

An important fact here is that the quark and gluon momentum fractions, i.e. (x)® and 
(x) 9 are empirically known with fairly good precision. For instance, the two popular PDF 
fits, i.e. MRST2004 and CTEQ5, give almost the same answer for (x)® and (x) 9 at least 
within the energy range Q 2 < lOGeV 2 . There also exist phenomenological fits for the 
longitudinally polarized PDFs, which contains the information on AE^ and Ag, although 
with larger uncertainties compared with the unpolarized case. These phenomenological 
PDFs can therefore be used for estimating the orbital angular momenta carries by the 
quark fields and the gluon fields through the relations : 

L Q = J Q -^AT, Q , (190) 
L 9 = J g -Ag. (191) 

The values of and L 9 at Q 2 = 4 GeV 2 estimated in this way are shown in Table IV1 Here, 
we use MRST2004 PDF fit to estimate J Q and J 9 0. On the other hand, A£ Q and Ag 
are estimated by using three independent PDF fits, i.e. LSS2005, DNS2005, and GRSV2000 
Here, all of the three independent fits we are using correspond to the MS 
regularization scheme. As one sees, there are sizable uncertainties for the phenomenological 
values of AS^ and Ag. Still, a common conclusion obtained from all these PDF fits is 
a very important role of the quark orbital angular momentum. The table |3 shows, at the 
least, that the magnitude of the quark orbital angular momentum is comparable with that 
of the intrinsic quark spin even at the scale of Q 2 = 4 GeV 2 . Since the quark orbital angular 
momentum is a rapidly decreasing function of the energy scale, while the scale dependence 
of AS^ is very weak, this means that the former is dominant over the latter at the scale 
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below Q 2 ~ 1 GeV 2 where any low energy models are supposed to hold. Naturally, the 
whole argument here is crucially dependent on one theoretical postulate that -6^,(0) ~ 0. 
Although it is supported by both of the LHPC and QCDSF lattice simulations, efforts to 
improve the accuracy of the lattice prediction should be continued, in consideration of its 
extremely important role in determining the quark-gluon contents of the nucleon spin. Also 
highly desirable is an analytical proof of it within the framework of (nonperturbative) QCD. 

TABLE V: The model independent predictions of the spin contents of the nucleon at Q 2 = 4 GeV 2 , 
based only upon one theoretical postulate -6^,(0) = 0. Here, all of the three independent fits for 
the longitudinally polarized PDFs correspond to the MS scheme. 



MRST2004 


LSS2005 


MRST2004 + LSS2005 


J Q J9 


AS Ag 




0.289 0.211 


0.198 0.368 


0.190 -0.157 




MRST2004 


DNS2005 


MRST2004 + DNS2005 


J® ja 


AS Ag 




0.289 0.211 


0.313 0.477 


0.133 -0.266 




MRST2004 


GRSV2000 


MRST2004 + GRSV2000 


J Q J9 


AS Ag 




0.289 0.211 


0.137 0.623 


0.221 -0.412 



V. CONCLUDING REMARKS 

In this paper, we have investigated the generalized form factors of the nucleon, which 
will be extracted through near-future measurements of the generalized parton distribution 
functions, within the framework of the CQSM. A particular emphasis is put on the pion 
mass dependence as well as the scale dependence of the model predictions, which we com- 
pare with the corresponding predictions of the lattice QCD by the LHPC and the QCDSF 
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collaborations carried out in the heavy pion regime around m v ~ (700 ~ 900) MeV. The 
generalized form factors contain the ordinary electromagnetic form factors of the nucleon 
such as the Dirac and Pauli form factors of the proton and the neutron. We have shown that 
the CQSM with good chiral symmetry reproduces well the general behaviors of the observed 
electromagnetic form factors, while the lattice simulations by the above two groups have 
a tendency to underestimate the electromagnetic sizes of the nucleon. Undoubtedly, this 
cannot be unrelated to the fact that the above two lattice simulations were performed with 
unrealistically heavy pion mass. 

We have also tried to figure out the underlying spin contents of the nucleon through the 
analysis of the gravitoelectric and gravitomagnetic form factors of the nucleon, by taking care 
of the pion mass despondencies as well as of the scale dependencies of the relevant quantities. 
After taking account of the scale dependencies by means of the QCD evolution equations at 
the NLO in the WS scheme, the CQSM predicts, at Q 2 = 4 GeV 2 , that 2 J Q ~ 0.68, A£ Q ~ 
0.32, and 2 L Q ~ 0.36, which means that the quark orbital angular momentum carries sizable 
amount of total nucleon spin even at such a relatively high energy scale. It contradicts the 
conclusion of the LHPC and QCDSF collaborations indicating that the total orbital angular 
momentum of quarks is very small or consistent with zero. It should be recognized, however, 
that the prediction of the CQSM for the total quark angular momentum is not extremely 
far from the corresponding lattice prediction 2 J® ~ 0.6 at the same renormalization scale. 
The cause of discrepancy can therefore be traced back to the LHPC and QCDSF lattice 
QCD predictions for the quark spin fraction AX^ around 0.6, which contradicts not only 
the prediction of the CQSM but also the EMC observation. As was shown in our recent 
paper A£^ is such a quantity that is extremely sensitive to the variation of the pion 
mass, especially in the region close to the chiral limit. More serious lattice QCD studies on 
the m^-dependence of AS^ is highly desirable. 

Worthy of special mention is the fact that, once we accept a theoretical postulate i?2o(0) = 
0, i.e, the absence of the net quark contribution to the anomalous gravitomagnetic moment 
of the nucleon, which is supported by both of the LHPC and QCDSF lattice simulations, 
we are necessarily led to a surprisingly simple relations, J® — ^ (x)® an d J 9 = \ (x) 9 , i.e. 
the proportionality of the linear and angular momentum fractions carried by the quarks 
and the gluons. Using these relations, together with the existing empirical information 
for the unpolarized and the longitudinally polarized PDFs, we can give model-independent 
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predictions for the quark and gluon contents of the nucleon spin. For instance, with combined 



use of the MRST2004 fit |64j and the DNS2005 fit |68|, we obtain 2 J Q ~ 0.58, AS Q ~ 0.31, 
and 2lfi ~ 0.27 at Q 2 = 4GeV 2 . Since (as well as J®) is a rapidly decreasing function 
of the energy scale, while the scale dependence of AS^ is very weak, we must conclude that 
the former is even more dominant over the latter at the scale below Q 2 ~ 1 GeV 2 where any 
low energy models are supposed to hold. 

The situation is a little more complicated in the flavor-nonsinglet (or isovector) channel, 
because B^ d (0) = G^(0)-G%$(0) = 2 J u - d - (x) u - d ^ 0, and also because the CQSM 
and the lattice QCD give fairly different predictions for (0). As compared with the 

lattice prediction for G^f^di^) around 0.8, the predictions of the CQSM turns out to be 
around 0.4. We have argued that the relatively small value of G^f^ (0) obtained in the 
CQSM is intimately connected with the small x enhancement of the generalized parton 
distribution E^~ 1 \x, 0, 0), which is dominated by the clouds of pionic qq excitation around 
x ~ 0. (We recall that the 2nd moment of (x, 0, 0) gives G^f^(0).) Unfortunately, 

such a ^-dependent distribution as E^~ l \x, 0, 0) cannot be accessed within the framework 
of lattice QCD. Still, the predicted small x behavior of E^J d (x 1 0, 0) = E^f l \x 1 0, 0) as well 
as of f u ~ d (x) indicates again the importance of chiral dynamics in the nucleon structure 
function physics, which has not been fully accounted for in the lattice QCD simulation at 
the present level. 
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APPENDIX A: PROOF OF THE MOMENTUM SUM RULE 

Here we closely follow the proof of the momentum sum rule given in |31|, by taking into 
account a necessary modification in the case of m n ^ 0. The starting point is the following 
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expression for the soliton mass (or the static soliton energy) : 

M N = N c S P [6(E -H + ie)H] - (H^H ) + E M , (Al) 

with 

E M [F(r)] = -fiml J [cos F(r) - 1] d 3 x. (A2) 

The soliton mass must be stationary with respect to an arbitrary variation of the chiral field 
U or equivalently the soliton profile F(r), which lead to a saddle point equation : 

Sv[9(E -H + ie)5H] + 5E m = 0. (A3) 

Here we consider a particular (dilatational) variation of chiral field 

U(x)^U((l+0x). (A4) 

For infinitesimal £, we have 

6U = U((l+£)x)-U(x) ~ £x k d k U(x), (A5) 

so that 

5H = M^°£x k d k U^ = £ [x k d k} M^IP 6 ] = £ ([x k d k , H] - r/V^fc)- (A6) 
Noting the identity 

Sv(0(E o -H + ie)[x k d k ,H}) = Sp([H,9(E -H + ie)]x k d k ) = 0, (A7) 

we therefore obtain a key identity 

£Sp[9(Eo-H + ie)(-i) 1 ° 1 k d k ] = -5E M . (A8) 

Now, by using (|A2|) together with the relations, 

SF(r) = £rF'(r), (A9) 
5 cos F(r) = - sin F(r)5F = -£r sin F(r) F'(r), (A10) 

we get 

rod 

8 Em = i-^flml \ drr 3 sin F(r)F\r) (All) 
Jo 

fOO d 

= -i-\-Kflm\ J drr 3 — cos F(r). (A12) 
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Here, taking account of the boundary condition 



F(Q) = n, F(oo) = 0, (A13) 



we can manipulate as 



drr —cos Fir) = / dr r d — (cos F(r) - 1) (A14) 
o dr Jo dr 

= r 3 (cosF(r) - - 3 / drr 2 (cos F(r) - 1) (A15) 

Jo 



3C 



= -3 / rfrr 2 (cos F(r) - 1). (A16) 
Jo 

We thus find an important relation : 

5£ M = -3££ M - (A17) 
Putting this relation into (|A8|) . we have 

£Sp(0(£ o + p) = 3££ M , (A18) 

or 

i Sp(#(£ -H + ie)a ■ p) = E M . (A19) 

If we evaluate the trace sum above by using the eigenstates of the static Dirac Hamiltonian 
H as a complete set of basis, (jA19|) can also be written as 

Y / (n\\a-p\n) = E M , (A20) 

which is the relation quoted in ([84)1 . We point out that our result has a correct chiral limit, 
since Em — » as — > and therefore 

lim V (n| - a • p In) = 0, (A21) 



in conformity with the proof given in ref 
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